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We study ring-exchange models for bosons or XY-spins on the triangular lattice. A four-spin 
exchange leads to a manifold of ground states with gapless excitations and critical power-law corre- 
lations. With a nearest-neighbour exchange, fluctuations select a four-fold ferrimagnetically ordered 
ground state with a small spin/superfluid stiffness which breaks the global Z7(l) and translational 
symmetry. We explore consequences for phase transitions at finite temperature and in an in-plane 
magnetic field. 
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I. INTRODUCTION 

Multi-spin exchange models incorporating ring- 
exchange processes have been of interest since the early 
studies of magnetism in solid Helium-3i*£. Ring-exchange 
processes could also play an important role in Wigner 
crystals near the melting densitjii^, and in Mott in- 
sulators which retain a fair degree of local charge 
fluctuations 5 . Indeed, neutron scattering experiments 6 
in insulating La2Cu04 have shown that aspects of the 
spin wave dispersion in the antiferromagnet may be un- 
derstood by invoking ring-exchange terms. Another rea- 
son for the interest in such models is that they may sup- 
port spin-liquid phases which are translationally invari- 
ant Mott insulators with no magnetic order, as indicated 
from numerics on triangular and Kagome lattices^. This 
has been established analytically in some models with 
U(l) symmetry 8 '-, which may also be viewed as boson 
models. Finally, many models such as the quantum dimer 
model on a triangular lattice^, the easy-axis version of a 
generalized Heisenberg model on the Kagome lattice?- and 
the easy-axis antiferromagnet on the pyrochlorc lattice 11 
may be mapped onto effective ring-exchange models in 
their low energy subspace of states. Thus, understand- 
ing the phases and phase transitions in ring exchange 
models is important. 

In this paper, we will focus on an XY ("easy-plane") 
ring-exchange model on the triangular lattice with both 
4-spin and 2-spin exchange terms (with strengths K and 
J respectively) in the regime J-Cif. This model is inter- 
esting from several points of view. First, several systems 
such as Wigner crystals and certain organic Mott insula- 
tors form a triangular lattice of spins with possibly appre- 
ciable ring-exchange processes, and the XY spin problem 
is one tractable limit of such models. Second, easy-plane 
magnets on triangular lattices are known to exist (though 
mostly as stacked layers forming a three dimensional sys- 
tem) and this could be of some relevance to them - 
in fact, these systems have motivated several studies on 
easy-plane triangular quantum Heisenberg models^. Fi- 
nally, earlier work by some of us has shown that re- 
lated models on the Kagome lattice are fractionalized 8 , 
whereas they support a critical phase (the "exciton Bosc 



liquid") on the square lattice^. This was essentially 
shown by perturbing around the pure ring-exchange limit 
(J = 0) where the Kagome and square lattices respec- 
tively have 0(L 2 ) and O(L) conserved quantities. It is 
thus worthwhile to examine cases where such extensive 
symmetries are absent from the outset even in the pure 
ring-exchange model, as is the case on the triangular lat- 
tice. 

The principal results of this paper are: (i) We show 
that the U(l) four-spin exchange model on the triangu- 
lar lattice (with S z — or half-filling for the bosons) 
has a manifold of degenerate ground states and the cor- 
relation functions are critical in the ground states. We 
identify the appropriate symmetries and conserved quan- 
tities which give rise to this, (ii) Perturbing in the 
nearest neighbor exchange (boson hopping) we find that 
fluctuations select a four-fold set of states from the de- 
generate manifold and the system develops long-range 
order at zero temperature. The ordered states break 
the global U(l) symmetry as well as translational sym- 
metry, (iii) In the ordered ground states, the super- 
fluid/spin stiffness ~ J 2 / K and is very small for small 
J/K. At any small nonzero temperature, the U(l) sym- 
metry is restored and we are in a phase with power-law 
phase/spin order, which gives way to a disordered phase 
via a Berezinskii-Kosterlitz-Thouless 14 (BKT) transition 
for T > Tbkt ~ J 2 /K. However, the discrete symme- 
try of the broken translations is not destroyed until a 
higher temperature T c ~ J. We analyze phase transi- 
tions at finite temperature and with in-plane magnetic 
fields within a Landau theory, and discuss the phase di- 
agram and experimental consequences. 

The outline of the paper is as follows. In Section II, we 
will present the Hamiltonian. In Section III we analyze 
its symmetries and conserved quantities, and define the 
dual representation of the model which allows us to an- 
alyze the instabilities of the system towards charge and 
energy ordering. We present arguments and numerical 
results to show that the model with J = has criti- 
cal power-law correlations at T — 0. In Section IV, we 
perturbatively analyze effects on introducing the near- 
est neighbor exchange (boson hopping) and find a phase- 
ordered state. We also discuss the phase diagram within 
Landau theory as a function of temperature and in-plane 



2 



(§ 7 r Sj 




© '->' 

(a) (b) 

FIG. 1: (a) The three kinds of plaquettes on the triangular 
lattice. The ring-exchange process involves the spins located 
at the points 1 - 4 of each plaquette. The labelling is chosen to 
coincide with that appropriate for the definition dual variables 
in the dual model discussed in Section 11(C), the sites P form 
the sites of the dual Kagome lattice, (b) Rotating <p — > ip + n 
on the sites indicated by open circles (which form one of four 
possible triangular sublattices) changes the sign of the ring- 
exchange term K in the Hamiltonian. Also shown are the 
basis vectors a, b for the triangular lattice. 



magnetic-field. We close with a discussion and specu- 
lations for three-dimensional generalizations and SU(2) 
invariant models in Section V. 



II. MODEL 

We define the model in rotor variables as 

H = K cos(y?i - ip 2 + <P3 ~ <P4) 
p 

+ J^2cos(<p i -<p j ) + — ^2(n i -n) 2 (1) 

where ipi is the phase of the boson variable (tpi — and 
ifi = 2n are identified) and n, is the the canonically con- 
jugate boson number (respectively the angle and angular 
momentum of a U(l) rotor), satisfying the commutation 
relation Upi, nA = iSi.j. The terms ip% . . . if4 in the first 
term of H denote angles around 4-site plaquettes of the 
triangular lattice — there are three such kinds of plaque- 
ttes as shown in Fig. 1(a). The J-term denotes nearest- 
neighbor hopping of bosons while U denotes a repulsive 
interaction between bosons. With U/K — > oo, we can 
identify S z (r) = n r — 1/2, S (r) = exp(±iip r ) and this 
model reduces to a S = 1/2 XY quantum spin model 
with n = 1/2 corresponding to total S z = 0. 

With J = 0, changing ip r — > cf> r = ip r + n on the sites 
indicated by open circles in Fig. 1 (b) changes the sign of 
K (there are four such choices of a triangular sublattice 
on which to make this transformation, as is clear from 
the figure). Thus, for J = 0, the ground state energy is 
independent of the sign of K. For nonzero J however, 
the sign of K is important. 



For K < 0, the leading perturbative corrections on in- 
cluding a nonzero J are O(J), and the energy and phases 
depend on the sign of J — for J < we get a ferromag- 
netic phase, while J > leads to a V3 x Neel order, 
which are also the ground states for large \J/K\ (where 
we may ignore the 4-spin term) . Thus there are no phase 
transitions at any nonzero J and we only obtain the well- 
studied phases. 

For K > 0, as appropriate for say spin degrees of free- 
dom in an electronic Mott insulator, it turns out that 
the leading perturbative corrections to the free energy 
are 0(J 2 ). To this order, the free energy is indepen- 
dent of the sign of J and we obtain a ferrimagnetic phase 
which breaks the global U(l) invariance as well as four- 
fold translational symmetry. The physics of this phase 
and the phase transitions out of it will be the focus of 
this paper. For large \J\/K we of course recover the con- 
ventional phases mentioned above. 

To make an estimate of the coupling constants in one 
case, let us imagine starting from the Hubbard model for 
electrons (with nearest neighbor hopping) at half-filling 
on a triangular lattice and perturbing in t/U to derive 
an effective spin model with ring-exchange terms in the 
insulator. This takes the form 

ffspin = *X [(Si • S 2 )(S 3 • S 4 ) + (Sr • S 4 )(S 2 • S 3 ) 

□ 

- (Si-S 3 )(S 2 -S 4 )]+^Ji J S i -S J - (2) 

id 

where the first term involves all 4-site plaquettes on the 
triangular lattice. Adapting results from Ref£ to this 
case, we find J = U 2 /U - 28t 4 /U 3 , J' = J" = 4t 4 /U 3 , 
K = 80t 4 /U 3 . For U/t = 6, numerical results^ show 
that the model is in an insulating phase though still close 
to the metal-insulator boundary; in this case J/t = 0.53, 
J'/t = J" jt = 0.02, K/t = 0.37. Clearly, it seems 
that one can describe the spin degrees of freedom in 
the insulator by setting the further-neighbor couplings 
J' = J" = and retaining only nonzero J, K. Indeed, 
exact-diagonalization studies^ find a spin-liquid phase for 
a closely related model in qualitative agreement with 
the Monte Carlo results on the Hubbard modeli^. It 
thus seems profitable to understand the above model for 
J/K<ti 1 as a starting point to analyze this full problem. 
For U/K — > co, the Hamiltonian in Eq.QJis precisely the 
XY limit of the above model obtained by setting terms 
containing S z to zero. 

III. MODEL WITH J = 

A. Symmetries 

For J = 0, we can change the sign of K by shifting 
fr —> </>r = <Pr + ^ on any one of four sublattices of 
the triangular lattice indicated earlier (we will work with 
these 'rotated' variables 4> r and a ferromagnetic 4-spin 
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term for convenience). This means shifting ip r — > (f> r + 7r 
on any two such sublattices leaves the action invariant 
and is a symmetry operation. We can identify this with a 
conservation of total boson number modulo two on alter- 
nate rows of the triangular lattice. These rows can run 
in any of the three symmetry directions of the lattice, 
this corresponds to a total of four symmetry operations 
including the identity. 

The rotor Hamiltonian above describes bosons hopping 
on plaquettes of the triangular lattice. For J = 0, the dy- 
namics conserves the center of mass of the bosons, since 
the bosons hop equal distances in opposite directions on 
any plaquette. In other words the center of mass oper- 



Rc.m. ) H 



= 0, and is a 



ator R c . m . = rh r satisfies 

constant of motion. Thus, in deriving the imaginary-time 
path integral for the partition function in terms of <f>, we 
may insert factors of exp(jq-R cm .) with an arbitrary vec- 
tor q, since R c . m . commutes with H. This corresponds 
to 4> r — > <j) r + q • r being a symmetry of the phase-action. 

If we identify cf> = with <p = 2ir, or keep track of 
the fact that the charge is quantized the first symmetry 
operation discussed above is a special case of the second, 
and corresponds to choosing q = Qi (i — 1,2,3) where 
the momenta are at the center of the edges of the 
Brillouin zone, indicated in Fig. 4(a). The above sym- 
metries for J = imply we cannot have terms which 
depend on (V^>) in the phase action (in the path integral 
for the partition function). The lowest allowed gradi- 
ent term on which the action can depend is of the form 
(S7 2 4>). We explicitly confirm this in a spin-wave calcula- 
tion in Section III.C, where we show the collective mode 
disperses as w 2 (k) ~ |k| 4 for small |k|, and calculate its 
consequences for correlation functions. The presence of 
a nonzero J explicitly breaks these symmetries, in fact 
J = will turn out to be a critical point of the model. 



B. Plaquette duality and dual action 

For J = 0, it is useful to define a dual representation for 
the model which permits us to obtain spin- wave theory as 
a well-defined limit of the model and to analyze instabil- 
ities of the spin-wave phase towards charge and energy 
ordering. For this purpose, we use a plaquette duality 
transformation^ and work with dual variables Op and 
Np which reside at the centers of each of the 4-site pla- 
quettes of the triangular lattice — these sites (labelled by 
P) form a Kagome lattice with nearest neighbour spacing 
1/2 in units of the lattice spacing of the original triangu- 
lar lattice. Alternatively, we may label sites on the dual 
lattice as (r, a), viewing the Kagome lattice as a triangu- 
lar lattice three sites (labelled a — 1,2,3) per unit cell, 
and we will use this notation whenever convenient. Thus, 
there are three times as many sites (and hence degrees 
of freedom) on the dual lattice compared to the original, 
and this is reflected in a redundancy in the description 
which we will discuss shortly. Define 
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FIG. 2: (a) The original triangular lattice (dashed lines) and 
the dual Kagome lattice (bold lines). The sites '+' and '-' 
indicate the sites of the Kagome unit around the site r of 
the triangular lattice, and are referred to as 'hex' and 'star' 
respectively in the paper, (b) The "gauge" transformation on 
the dual lattice — shifting by e on the indicated sites as 
shown, with e/n being an integer, is an exact symmetry of 
the dual Hamiltonian. 
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where the angles (f>\ ... 4>A in Eq. @ are around plaquette 
P of the lattice with the convention indicated in Fig. 1(a). 
The sites hex and star lie on the 12-site unit of the 
Kagome lattice around the site r as shown in Fig. 2(a) 
(marked by '+' and '-' respectively), while the center of 
each Kagome unit lies on the triangular lattice (and is 
labelled r). In these variables, the dual Hamiltonian on 
the Kagome lattice is 



#dual = if ^COs(TTiVp) 



U 
2^2 



(5) 
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and the dual variables satisfy [iVp,#p/] = i8p t p', with 
dp I having integer eigenvalues, while Np has a contin- 
uous spectrum. 

We pause to elaborate further upon the redundancy 
of the variables. Since only the combination n r of the 
Op variables is physical, the Hamiltonian will be invariant 
under any integer multiple of 7r shift of the 9p that leaves 
all the n r invariant. This can be done locally as follows. 
Pick any hexagon of the dual Kagome lattice, and choose 
two neighboring sites on this hexagon. Shift Op — ► Op + ir 
on those sites, and simultaneously take Op — ► Op — it on 
the sites on the opposite side of this hexagon, leaving the 
remaining two site of the hexagon (and all the other sites 
not on this hexagon) untouched as shown in Fig. 2(b). 
Since this leaves all physical properties invariant, this 
should be regarded as a local gauge symmetry. Such 
transformations are generated by the unitary operators 



Gr(K Q ) 



ac, 



(6) 



where K Q with a = 1,2,3 are the three reciprocal lattice 
vectors for the triangular lattice. In the spirit of a gauge 
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theory, we could restrict our consideration to physical 
gauge- invariant states which satisfy G r (K a ) = 1. This 
can be rewritten suggestively as 

J2( r P -r)N P = O(moda), (7) 

hex 

where a is an arbitrary Bravais lattice vector for the tri- 
angular lattice. As in Refs£ji£, the operator Np may 
be regarded as the local "vortex number" (modulo 2) 
on the site P of the dual lattice. Thus the condition of 
Eq. J7J requires that the local vortex center of mass on 
each hexagonal plaquette of the Kagome lattice vanishes 
- the ambiguity by a Bravais lattice vector corresponding 
exactly to the ambiguity in Np by a shift of 2. Thus the 
gauge-invariance of the dual description is related to the 
immobility of "vortices" in this model. 

Using a Trotter decomposition, we may write the par- 
tition function as a discretized imaginary-time path in- 
tegral in the standard manner. This leads to Z = 
f V9(t,t) exp(-S [9(r,r)\) with the action 

^ 2 

Sduai = — 5-hi( — T-) y (9 raT — 9 r aT+ i) 
ir z e T K z — ' 

r,a,r 

+ gE(E^-E^--) 2 (8) 

r hex star 

where the field 9 rta-T /n is an integer- valued field, e T = 
j3/N T with /3 = 1/T is the inverse temperature, and N T 
is the number of imaginary time slices. The quantum 
problem at a temperature T corresponds to the contin- 
uum limit e T — > 0, N T — ► oo with fixed l/e T N r = T. 
Thus, we reduce the quantum model to an anisotropic 
(2 + l)-dimensional classical model — since there is no 
sign-problem, this proves useful for carrying out Monte 
Carlo simulations in the dual representation to numer- 
ically test for charge ordering and plaquette-energy or- 
dering in the ground state at arbitrary U/K. 

C. Effective theory and spin-wave approximation 

We may rewrite the dual action as 

Sdual = o hl( — ) ^ {9 r ,a,r — #r,a,r+l) 2 

r,a,r 

r hex star 

- E V% q COs(2q9 r<a , T ) (9) 
<j,r,T,a 

where the bare couplings v\ q are chosen to enforce the 
integer constraint on 9, and may be singular. We will 
now write down a low energy effective action in terms of 
$r,a,r = W\f — (Qmn)/2. Here [9]f symbolically denotes 
an integration over the fast (high-energy) degrees of free- 
dom while retaining the spatial structure of the lattice, 



Q = 2-Kn, and the site r = ma + nb in terms of the basis 
vectors a, b of the triangular lattice shown in Fig. 1(b). 
The factor we have subtracted out eliminates the mean 
density n from the second (boson repulsion) term in the 
Hamiltonian. 

If we ignore the terms with v% q , we obtain a Gaussian 
action in terms of $. Studying this action motivates us 
to guess that the effective low energy action may be of 
the form Sf^=sfl x + 5« al with 

(10) 

and 

^dual = / E V2 1 cos ( 2 <?^r,Q,r + Qqmn) + Shighcr (11) 

where we have defined renormalized couplings £/(k), /C(k) 
(which are smooth nonzero functions of k) and v^q, the 
site r = (to, n) as before, and the 3x3 matrix G a ^(k) = 
A* a (k)A!3{k) with 

Ai(k) = 1 + e- ik « - e - lk «- lkb - e ik " 

A 2 (k) = 1 + e - lk «- lkb - e~ lka - e- lkb 

A 3 (k) = 1 + e-' lkb - e lka - e -^- 2fe o, (12) 

where k a = k • a and fc& = k • b. S'higher denotes other 
terms which might be generated in deriving the low en- 
ergy action. For the present purposes, we will not write 
down the explicit form of these terms but we will discuss 
them later. 

If we set V2 q = and ignore S'higher, we can obtain the 

eigenmodes of the Gaussian action Sj^L by diagonalizing 
G a p{k). We find there is one nonzero eigenvalue A(k), 
and two zero eigenvalues. This corresponds to one dis- 
persing mode (Si) and two nondispersive modes (62,3) 
with zero energy (flat bands) and the action at this level 
takes the form 

^ = /E(^)ie,<^ 

+ /" k .S. 3 (^w) |eo(k " )|2 ' (13) 

where we have set A(k) = £ 2 (k)/W(k)/C(k) so that £(k) 
appears as the excitation energy of the collective mode 
0i. Explicitly, 

A(k) = 2 [6 — 3cosfc a — Scoskb — 3cos(fc a + kb) (14) 
+ cos(fc a - kb) + cos(2fc a + kb) + cos(2fc b + k a )], 

thus the dispersing mode ©i has one gapless point at 
k = (0,0). The leading dispersion away from this point 
is £ 2 (k 0) = (9/16)W(0)/C(0)|k| 4 — as stated in sub- 
section A, this dispersion is dictated by the symmetry of 
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conservation of center of mass of the bosons which forbids 
a term ~ |k| 2 . 

We can also now see that the extra degrees of freedom 
(two per site of the triangular lattice) introduced in go- 
ing over to the dual description manifest themselves as 
zero energy modes corresponding to unphysical fluctua- 
tions, while the dispersing mode corresponds to physi- 
cal fluctuations. In fact, by Fourier transforming a gen- 
eral superposition / 2 (k)02(k) + /3(k)9 3 (k), we can show 
that the zero modes indeed arises from the local gauge 
symmetry generated by Eq. (j^J in real space. More pre- 
cisely, in the coarse-grained quadratic low energy theory, 
these discrete symmetry operations are promoted to con- 
tinuous shifts (with arbitrary |K Q |). Similarly, Fourier 
transforming / 1 (k)0 1 (k), we find that the physical fluc- 
tuations are composed precisely of the gauge invariant 
boson densities. 

At this stage, we may reintroduce the boson phase 
fields by a Hubbard Stratonovitch decoupling of the time 
derivative term of the physical fluctuation ©i (as in 
Rcf. 13 ), namely 

e -*£®\e>Mf = J w . e -£$l*M a (15) 

x e ^lb) (<A(k)e?(k) -^ (k)ei(k)) , 

Provided V2 q and Shigher are irrelevant and can be ig- 
nored, we can integrating out the field ©i which leads to 
a Gaussian action for the boson field <f>, 

This corresponds to a spin-wave (harmonic) approxima- 
tion and may be used to evaluate boson correlation func- 
tions to zeroth order in Vi q . 

D. Boson correlation functions in spin- wave theory 

To characterize the Gaussian theory , we evaluate 
the space-time correlation functions for the boson cre- 
ation operator exp(i</> rT ). At long times or large separa- 
tions, these correlation functions reduce to 

i / »([>) 

/g'0r,Og-«>D,O\ ^ | r | 2*-v/3 V K <°) 
1 / W(0) 

( e <*O,r e -i*O.0) „ | T | wsV*^. (17) 

Clearly, the Gaussian theory describes a critical liquid of 
the bosons with power law correlations arising from the 
gapless excitations dispersing as £ (k) ~ |k| 2 . From the 
long time behavior of the two-point correlation in a finite 
size system, we may evaluate the finite size scaling of the 
gap to adding a particle to be just A(L) ~ U(0)/L 2 . 
In the thermodynamic system, the low energy density of 



states per unit volume of collective excitations is given 
by 

N(u) = ~£*( W -£(k)) 

k 

which is a constant depending on the interactions. Thus, 
all low-energy long-wavelength properties of the liquid 
are determined in terms of the k — > behavior of the 
functions i/(k),/C(k). We next present arguments and 
numerical evidence that this Gaussian description of a 
critical liquid in the model with J = may be valid even 
in the presence of terms V2 q and Shighcr- 



E. Argument for irrelevance of higher order terms 

To see whether the term v^q is relevant, we can evaluate 
the correlation functions of this operator, and we find it 
is local in space and exponentially decaying in time, 

(cos2$ riQir cos2tV,3,o) ~ 5 rfi 5 a ^ exp(-7|r|) (19) 

with a constant 7 > 0. This is not hard to understand — 
it arises from the fact that the variable 1? is a combination 
of the eigenmodes 61,2,3 and cos2i9 r , QiT is not "gauge- 
invariant" . As a result its spatial correlations are local, 
while the fluctuations of the zero-modes 62,3 determine 
the exponential temporal decay. Clearly, vi q is irrelevant. 

Consider the possible forms we can obtain for Shigher- 
These terms would be of the form of cosines involving 
multiple $'s at different space-time points. However, sim- 
ilar arguments as above apply to these operators. The 
only operators which can have nonzero correlations at 
nonzero separation would be gauge-invariant combina- 
tions of the i?'s. As we have shown, these are the local 
boson densities. Thus, if we admit only slow fluctua- 
tions of 0i, these would take the form of weak density- 
density interactions within a perturbative treatment and 
could renormalize the coefficients of the Gaussian action, 
but not cause an instability. Of course, such density- 
density interactions might lead to a charge-ordered state 
at strong couplings, but our arguments show that the 
Gaussian theory is perturbatively stable^. 



F. Numerical results 

To confirm the above arguments we have carried out 
Monte Carlo calculations on the dual model using a 
Metropolis algorithm, for lattice sizes upto 432 spatial- 
sites and 48 time-slices, and periodic boundary condi- 
tions. We find no evidence of any charge ordering for 
h = 1/2, even for large U/K. We have compared 
the density-density correlations and find agreement with 
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FIG. 3: The dependence of the scaled compressibility (e T K)n 
on the bare boson repulsion U, for various e T K. Within spin- 
wave theory, we expect (e T K)K, = K/U(Q) as discussed in the 
text. From the figure we see that K/U(0) coincides with K/U 
for small U (as indicated by the dashed line), but deviates 
and tends to a constant K/U c s ~ 0.2 as the bare U — » oo (as 
shown by the solid line drawn as a guide to the eye). Typical 
error bars are indicated at one point. We thus expect the 
S = 1/2 XY spin model with J = to be compressible, with 
W(0) = U cS . 



what we expect from a Gaussian theory «5^ ual . In par- 
ticular, we show results for the compressibility k defined 
through k — Xrm(k — > 0, 0J n — 0) where 

Xnn(k,o,„) =^ e - 4k - r +-" T (n(r,r)n(0,r)). (20) 

r.r 

Within the Gaussian theory, k = l/e T U(0). In Fig. 3 we 
plot the scaled compressibility (e T K)n which is expected 
to scale as K/U(0). We find that that it is nonzero 
at all U/K in the quantum (r-continuum) limit, with 
U(0)/K ~ 5 as U/K — > oo. We have also checked 
and found no evidence for energy ordering on the pla- 
quettes of the triangular lattice. Thus, we believe that 
for n = 1/2, the ring-exchange model with J = is well 
described by the Gaussian fixed-point action S^uLi where 
terms with V2 q and all other terms Shighcr are irrelevant. 



IV. PERTURBING IN NEAREST NEIGHBOR 
COUPLING J 

A. Ordered ground state 

To perturb in the nearest-neighbor exchange J, it is 
convenient to work in the phase representation which we 
have argued remains a valid 'fixed-point' description of 



the pure ring-exchange model for all values of the bare 
coupling U/K. In this, we make the approximation of 
completely ignoring irrelevant operators — this seems 
reasonable in retrospect since we will find that the system 
acquires long-range phase order at T = 0. 

The nearest-neighbor perturbation we add to the 
Gaussian phase action takes the form 

S,J = jy^ ldT [(-1)" COs(0 m ,„ iT - 4> m +l,n,r) 
m,n 

+ ("I)™ COs(0 m ,„, T - <j>m,n+l,r) 

+ {-l) m+n COs(<p m , n:T - <p m+1 , n+hT )] (21) 

where we have chosen to label sites on the triangular 
lattice as r = ma + n&, and we are working in rotated 
variables 4> T for which the ring-exchange term is ferro- 
magnetic. 

To carry out a systematic expansion, we define 4> r — 
4> r + Q ■ r with an arbitrary vector Q. This allows us 
to examine at the effect of perturbing in J in any of 
the ground states of the degenerate manifold to see if 
fluctuations select any state. Perturbing in Sj around 
the pure ring-exchange Gaussian theory, we find that the 
corrections to the free energy vanish at leading order in 
J. At C(J 2 ), we find a correction to the free energy 
density 

Sf = -\PJ 2 (cos 2 Q a + cos 2 Q b + cos 2 (Q a + Q b )) 

X *^2 ( _1 )"( COS (^0,0,0 - ^1,0,0 + 4>7n,n : T + </>m+l,n,r))o 
m,n ^ T 

(22) 

where (. ..)o is the expectation value evaluated in the 
Gaussian theory, Q a = Q • a, Qb = Q ■ b, and we have 
used the six-fold symmetry of the triangular lattice to 
simplify certain intermediate expressions. The expecta- 
tion value may be analytically evaluated at small U / K 
by expanding the cosine term since its argument is small 
in the Gaussian theory, and we find that the sum over 
to, n is positive. Thus, (i) the corrections to the free en- 
ergy are independent of the sign of J, and (ii) minimizing 
6f with respect to variations of Q corresponds to maxi- 
mizing (cos 2 Q a + cos 2 Qb + cos 2 (Q a + Qb)) which leads 
to ordering at four wavevectors Q Q (a — 1 ... 4) which 
are shown in Fig. 4(a). We expect this order to per- 
sist even at large effective (U/K) (recall that the bare 
U/K — > oo in the 5=1/2 spin limit). To obtain the 
ordered states in the original phase variables, we rotate 
back 4>(r) ip(r) = <f>(r) — tt on the triangular sublattice 
and we then find, as shown in Fig. 4(b), that the result- 
ing ordered states break (i) global U(l) spin rotational 
invariance and (ii) four-fold translational invariance. For 
small J/K, such states have been shown to occur at fi- 
nite magnetic fields for the classical version of a related 
SU(2) invariant ring-exchange modeUi, and this ordering 
has been called uuud (for 'up- up- up- down' which is the 
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spin order on a 4-site plaquette, but in the plane in our 
case). Here, we find such an ordered state is stabilized 
by quantum fluctuations in the XY limit even in zero 
magnetic field. 



B. Fluctuations in the ordered state 

As shown above, in the presence of a nonzero J, the 
model has an ordered ground state. This order breaks 
U(l) invariance. At the same time, it also breaks trans- 
lational symmetry. In the boson language, this state may 
be viewed as a superfluid coexisting with broken transla- 
tional order. This broken discrete symmetry is also evi- 
dent from our implicit choice of a state Q Q to minimize 
the free energy. Thus, we expect two kinds of excita- 
tions at low energies about the ordered state — (i) phase 
fluctuations, and (ii) domain walls in the discrete order. 

It is easy to see that with <j) r — Q • r + <j> r , there is 
an energy cost ~ (Q — Q Q ) 2 for small deviations |(Q — 
Qq)| *C 1 where Q Q denotes any one of the ordering 
vectors. This implies that the long wavelength effective 
theory of phase fluctations around each of these ordered 
states is described by an action of the form 



(. 



1 



S = J ^—(dr<p) 2 + J s (Vv) 2 +K cS {V 2 <p) 2 j (23) 

with a nonzero J s , and effective couplings U e e and Aeff- 
For small U/K, J/K, we find U c s ~ U, K e g ~ K and 
J s ~ J 2 /KyV / K. Thus, phase fluctuations are con- 
trolled by the phase stiffness J s . 

A sharp domain wall separating regions with different 
Q Q order would cost an energy ~ K per unit length. 
However, since the model is spin-disordered for J = 0, 
we expect that a smooth deformation of the spin con- 
figuration by making a domain wall with nonzero width 
£d, would cost less energy. Indeed, using the above ac- 
tion one can estimate the energy per unit length of a 
straight domain wall as edom = JiKetf/^D + 72</s£d. 
where 71.2 ~ 1 are constants. Minimizing the energy cost 
with respect to £rj, we find that ~ W K e g/J s , and thus 
domain wall excitations cost an energy e*j om ~ \/K c gJ s 
per unit length. For small U/K, e* dom ~ J (U/K) 1 / 4 . 



C. Phase transitions at finite temperature 

At high temperature, we expect the model to be in 
a fully disordered phase where the U(Y) symmetry as 
well as the discrete broken symmetry is restored. To 
study these phase transitions, we appeal to the above 
estimates of the spin stiffness and domain wall energies. 
For U/K ~ 1, since the spin stiffness J s ~ J 2 /K, while 
the domain wall energy per unit length edom J 3* J s , 
we expect the U(l) symmetry to be restored via a BKT 
transition^ at a temperature Tbkt ~ Js once vortices 
are included in the effective theory, while the discrete 




^0 



(a) 




(b) 



FIG. 4: (a) The four ordering wavevectors Q Q (a = 1 ... 4) 
obtained by perturbing in the nearest neighbor coupling J, 
with the hexagon indicating th first Brillouin zone of the tri- 
angular lattice, (b) The phase/spin order in the original vari- 
ables in one of the four broken symmetry ground states. 
The order may be most easily thought of as alternating rows of 
ferromagnetic and antiferromagnetic spins. The bold dashed 
lines indicate bonds on which the spins point in the same 
direction, while the other bonds have antiferromagnetic spin 
orientation. As discussed in the text, these correspond to en- 
ergy ordering in the ground state with an associated four-fold 
broken discrete symmetry. 



symmetry is not restored until a much higher tempera- 
ture T c ~ J. Thus there is an intermediate phase with 
exponentially decaying spin correlations, but with the 
discrete order still present. This bears resemblance to 
proposed scenarios of chiral ordering in the XY antifer- 
romagnet on a triangular lattice as well as the SU(2) 
multiple-spin exchange model. However, in contrast to 
the XY model where a single scale, namely the two-spin 
exchange coupling, sets the scale for both transition tem- 
peratures, the transition temperatures are governed by 
distinct energy scales in the present case, and the dis- 
crete symmetry corresponds to broken translational in- 
variance. 

In order to better understand the discrete transition, 
let us identify the appropriate order parameter and con- 
struct a Landau theory. Since the discrete symmetry 
is easily identified with translations in the original vari- 
ables tp r , we will discuss it in terms of that. As shown in 
Fig. 4(b), the ordered state at T = has broken transla- 
tional invariance associated with ordering of cos(<fi — tpj ) 
on nearest-neighbor bonds. Since there is a nearest- 
neighbor spin exchange term J, this leads to ordering in 
the energy density. We identify the broken discrete sym- 
metry with energy ordering on the bonds. To identify the 
appropriate order parameter for this transition, note that 
we may write K r i — (S r • S r+ ai) where hi (i — 1, 2, 3) are 
the unit lattice vectors making an angle 27r/3 with each 
other (specifically, h\ — a, a-i — b, and S3 = — (a + b)). 
Labelling the four ordered states by n = ... 4 and with 
ao = 0, we can write the expectation value in any of these 
ordered states as — cos(Q, • r)cos(Qi • a M ). For a 
general superposition, 



E 

^=0...4 



cos(Qj • r) cos(Q l • a M ) (24) 
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and we identify i? M as the appropriate order parameter 
for the Landau theory. Under shifts of — ► + A, 
the physical correlation K r i is unchanged. Thus we re- 
quire the Landau functional to be invariant under such 
shifts. Further, studying the transformation of K T i un- 
der symmetry operations of the lattice (unit translations, 
7r/3-rotations, and reflections about the 3 mirror planes), 
we find that such symmetry operations correspond to all 
possible permutations of the set {B^}. For the Landau 
theory we find beyond quadratic order, one cubic invari- 
ant and two quartic invariants. It appears likely that the 
finite temperature phase transition restoring this discrete 
symmetry is in the same universality class as a 4-state 
Potts model. 

To summarize, we expect the system to exhibit with 
increasing temperature a BKT transition with Tbkt ~ 
J 2 /K from a state with power law spin correlations to a 
state with exponentially decaying spin correlations. The 
discrete broken symmetry associated with energy order- 
ing on the bonds of the lattice is expected to be restored 
above a temperature T c ~ </, with the transition being 
in the universality class of the 4-state Potts model. In a 
weak applied in-plane magnetic field, the U(l) invariance 
is lost, and the BKT transition would become rounded 
but the discrete transition would survive as a finite tem- 
perature transition. At large enough fields, we expect a 
first order transition to a state where all the spins point 
in the same direction and the translation invariance is 
restored. 



V. DISCUSSION AND CONCLUSIONS 

Studying an XY model with four-spin couplings and 
small two-spin exchange, we have obtained a ferrimag- 
netically ordered ground state and analyzed phase tran- 



sitions out of this phase. Since the ground state is phase 
ordered, we expect our result to be stable to small out-of- 
plane couplings involving S z . Similarly, introducing weak 
couplings between such ordered two-dimensional planes 
would lead to a three-dimensionally ordered stated, 
where the spin/phase order would persist to finite tem- 
perature. We thus expect this state might be of some 
interest for Mott insulators on a triangular lattice. If 
such an ordered state exists, it could be detected in neu- 
tron diffraction studies. The spin-disordered state with 
persisting discrete broken symmetry may be indirectly 
observable through lattice distortions if the spins couple 
to the lattice. 

We can also analyze the opposite limit from the one 
studied in the paper, where we assume "easy axis" 
anisotropy in the model in Eq. @ and weak J/K. In 
this Ising limit, we find eight degenerate ground states 
for the spin model, again of the uuud type, with the 
spins aligned along S^-axis. These states correspond to 
charge-ordered incompressible states for the bosons at 
density fl = l/4orfz. = 3/4 (four ground states at either 
density). It may be that any anisotropy which takes us 
away from the SU(2) invariant model leads to ordered 
ground states, while the S'[/(2)-invariant model seems to 
be a uniform spin liquidii^. It would be interesting to 
further explore this possibility. 
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